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ABSTRACT 



Context, Since 1995, more than 500 extrasolar planets have been discovered orbiting very close to their parent star, where they 
experience strong tidal interactions. Their orbital evolution depends on the physical mechanisms that cause tidal dissipation, and these 
are still not well understood. 

Aims. The purpose of this paper is to refine the theory of the equilibrium tide in fluid bodies that are partly or entirely convective, to 
predict the dynamical evolution of the systems. In particular, we examine the validity of modeling the tidal dissipation by the quality 
factor Q, as is commonly done. We consider here the simplest case where the considered star or planet rotates uniformly, all spins are 
aligned, and the companion is reduced to a point-mass. 

Methods. We expand the tidal potential in Fourier series, and express the hydrodynamical equations in the reference frame, which 
(-H . rotates with the corresponding Fourier component. The results are cast in the form of a complex disturbing function, which may be 

implemented directly in the equations governing the dynamical evolution of the system. 

Results. The first manifestation of the tide is to distort the shape of the star or planet adiabatically along the line of centers. This gen- 
erates the divergence-free velocity field of the adiabatic equilibrium tide, which is stationary in the frame rotating with the considered 
Fourier component of the tidal potential; this large scale velocity field is decoupled from the dynamical tide. The tidal kinetic energy 
C/2 ' is dissipated into heat through turbulent friction, which is modeled here as an eddy-viscosity acting on the adiabatic tidal flow. This 

dissipation induces a second velocity field, the dissipative equilibrium tide, which is in quadrature with the exciting potential; it is 
responsible for the imaginary part of the disturbing function, which is implemented in the dynamical evolution equations, from which 
one derives characteristic evolution times. 

Conclusions. The rate at which the system evolves depends on the physical properties of tidal dissipation, and specifically on how the 
eddy viscosity varies with tidal frequency and on the thickness of the convective envelope for the fluid equilibrium tide. At low fre- 
quency, this tide retards by a constant time delay, whereas it lags by a constant angle when the tidal frequency exceeds the convective 
. turnover rate. 

C*~) • Key words. Gravitation - Hydrodynamics - Stars: binaries: close - Stars: planetary systems - Stars: rotation - Planets and 

■ satellites: general 
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1 . Introduction and general context have allowed an improved quantification of the tidal dissipation 

in Jupiter (Lainey et al. 2009) and Saturn (Lainey et al. 2012). 

Stellar binaries provide precious information on the structure of xhis rene wed interest for the gravitational interaction between 

V ; stars and the dynamical processes operating in them. Analysis various ce l es tial bodies motivates us to refine our knowledge of 

5^ , of the motion of a binary's components allows to determine the physical proC esses that intervene in tides. 

■ their masses and radii, and through the apsidal precession it 

puts constraints on their internal mass distribution. Furthermore, The dvnarnica i evolution of a binary system is driven by the 

their orbital evolution due to tidal dissipation sheds some light convers i on of its mechanical energy into heat. Provided the 

on their formation and past history. Quite naturally, similar system loses no angular momen tum, it tends to the state of 

benefits may be expected from the study of extrasolar planets. minimum en ergy in which the orbits are circular, the rotation 

Indeed, since 1995, more than five hundred of those have of the components is synchronized with the orbital motion, and 

been discovered orbiting very close to their host star, and their the spins ^ aligne d. However, in very close systems with high 

number keeps increasing thanks to dedicated space missions mass ratio> such final state cannot be achieved: instead; the 

(CoRoT and KEPLER) and vast ground-based programs, with planet spirals toward the star and may even tually be engulfed 

powerful instruments such as HARPS (for a review, see Udry & by it (Hut 1981; Levrard e t al. 2009). To predict the fate of a 

Santos 2007 and references therein). Furthermore, recent s tudies binary systerri; one has t0 identify the dissi p ative processes that 

on the orbital evolution of natural satellites in our solar system achieve me conversion of kinetic en ergy into thermal energy, 

from which one may then draw the characteristic times of 

Send offprint requests to: F. Remus circularization, synchronization and spin alignment. Before 
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reviewing these processes, let us recall the two types of tides 
which operate in stars and in the fluid parts of giant planets. 



The equilibrium tide designates the large-scale flow induced 
by the hydrostatic adjustment of the star in response to the 
gravitational force exerted by the companion (Zahn 1966; 
we shall refer to that paper as Z66). On the other hand, the 
dynamical tide corresponds to the eigenmodes (gravity, inertial, 
or gravito-inertial waves) that are excited by the tidal potential 
(Zahn 1975; Goldreich & Nicholson 1989; Ogilvie & Lin 
2004; Goodman & Lackner 2009; Rieutord & Valdetarro 2010; 
Barker & Ogilvie 2010). These tides experience two main 
dissipative mechanisms: turbulent friction in convective regions 
and thermal dissipation acting on the gravito-inertial modes 
excited in radiative zones. 




We shall focus here on the equilibrium tide acting in the con- 
vective envelopes of solar-like stars or giant planets. The phys- 
ical description of this tide, with turbulent dissipation, has been 
formulated by Z66. Several authors have applied this theory to 
binary stars (see Koch & Hrivnak 1981; Bouchet & Zahn 1989; 
Verbunt & Phinney 1995), and they found that it agreed well 
with the observations (pre-main-sequence binaries, red giants, 
etc.), except for main-sequence stars older than about one Gyr 
(Mathieu et al. 2004). Here, we revisit the theoretical formalism 
of Z66 by expressing the velocity field of the equilibrium tide 
in the appropriate reference frame, i.e. that which rotates with 
each Fourier component of the tidal potential. This treatment 
filters out properly the dynamical tide, suppressing the spuri- 
ous 'pseudo-resonances' encountered in Z66. Our second goal, 
in the present work, is to improve the description of the equi- 
librium tide by tying the quality factor Q, which is commonly 
used in planetology and celestial mechanics, with the physical 
process of convective turbulent dissipation (Goldreich & Soter 
1966, Goldreich & Nicholson 1977). We shall assume here, for 
simplicity, that all spins are aligned, i.e. that the rotation axes 
of star or planet are perpendicular to the orbital plane; the case 
with non-aligned spins will be treated in a forthcoming paper. 
Furthermore, the differential rotation, which is present in stellar 
and planetary interiors is not accounted for in this first work and 
the rotation is taken as uniform. 



2. Reference frame and tidal potential 

2. 1 . The reference frame 

We consider a system consisting of two bodies A and B, of mass 
niA and ot#, and we undertake to describe the tide exerted by B 
on A, which we assume to be in fluid state, i.e. a star or a gi- 
ant planet. Due to their mutual attraction, they move in elliptic 
orbits around their common center of mass, but it is often con- 
venient to choose a reference frame Ka whose origin is placed at 
the center of A (designated by the same letter) and whose axes 
{X, Y, Z) have fixed directions in space. We assume that all spins 
are aligned along the (AZ) axis; therefore the position of B is en- 
tirely determined by the three following keplerian elements : the 
semi-major axis a of the relative orbit of B around A, its eccen- 
tricity e, and its mean anomaly M - nt where n is the mean 
motion. 



Fig. 1. Spherical coordinates system attached to the reference 
frame %a centered on A. The spin of Ais perpendicular to the 
orbital plane of B. The dashed line illustrates the orbit of B. 



2.2. The tidal potential 

The tidal force exerted by B on A derives from a potential, which 
may be expanded in Legendre polynomials as follows : 



s>2 



(1) 



where = (AB, AM), M being the current point of position r, 
and d the distance AB. Kepler's third law Q (iiia + to#) — n 2 a 3 
allows us to rewrite ( 1) in a slightly different form: 



U E 



(,,,) = ^^V^P l(c o S 0) 

m A +m B d s+l 



(2) 



Finally, introducing the angular variables (9, tp) related to the 
reference frame %a an d the associated Legendre functions 
P™ (cos 9), we get 



u E (r,o,<p,t) 



m B 



rriA + m B 



(s — to)! 



(s + m) ! 

(2 - 6 , m ) K (0) K (cos 9) exp [imp] f^Y" exp [-imv] 



m=0 

VS+1 



(3) 



where v designates the true anomaly of B. The next step is to 
expand the time-dependent terms in Fourier series: 



(a\ s+l v-i 
Yd) ex P[ _imv ] = 2-i G s,^,l-m(e) exp[-i7nf] , (4) 



where we made use of Kaula's functions Gj tPi q(e) (Kaula 1962). 
Furthermore, we have: 



P'f (0) = 



, (to + s)! 

(-l)^ — — — , if (s + to) is even 

2 S (^)!(^)! ■ (5) 

, if [s + m) is odd 



We chose the Kaula's functions in the Fourier expansion, instead 
of the Hansen coefficients, to prepare the ground for a subse- 
quent paper where we shall treat the case of non-aligned spins. 
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These functions G SM (e) are polynomials in e whose lowest 
power is \q\. We give their expressions in Table 1 for the most 
frequently used triplets {s,p, q), since for small eccentricity it 
will suffices to retain only a few q values (q e [[-2, 2 J). 
From now on, we shall keep only the time-dependent part of the 
potential Ue, which we call U, the tidal potential: 

U (r, 6,<p,t) = Y Yj Yj U ™ (r) P ™ (cOS 0) COS ~ lnf) ' 

(6) 



s>2 m>l I 



where the radial parts U" 1 of the potential are given by 

(s — m)\ 



2 s (^)i as 



G s ™ i- m (e) ; 



the tidal potential is thus the interaction potential, from which we 
have subtracted the acceleration responsible for the movement of 
the center of mass of A: Qnipld [l + (r^ ■ r) /of 2 ]. 
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Table 1. Values of the eccentricity function G JiM (e) in the 
quadrupolar case where s — 2. Note their symmetry property 
G s ,p,q{e) = G SiS - p - q (e). (FromLambeck 1980) 



We may further simplify this expression of the tidal potential if 
we consider binary systems that are separated enough to allow 
the companion to be treated as a point mass; this approximation 
has been discussed by Mathis & Le Poncin-Lafitte (2009). In this 
quadrupolar approximation, we shall retain only the lowest term 
of the potential (s = 2). Furthermore, since we are considering 
the coplanar case, the sole value of m that still remains is m = 2 
in Eq. (6) and U thus becomes 

U = -\ — — — nV V Gooi- 2 (e)Pl(cosG)cos(2<p-lnt). (8) 
4 niA + mg z j J ' ' 

3. The tidal velocity field 

3.1. Governing equations 

The tidal potential (8) induces in the star (or the giant planet) 
pressure and density perturbations and a modified velocity field 



W, which obeys the following equations of motion: 
DW 



P— = -VP-pV(<D+{/) + F vlsc , (9) 



of continuity: 



Poisson: 



Dp 

— + pV • W = , 
Df H 



V 2 = AnQp , 



(10) 



(ID 



and of entropy: 

Cpe-VH 

Ui> (for a radiation zone) (12a) 

p dH t 

(cste (in a convective zone) (12b) 



(7) where we made use of the total derivative 



— = — + (W • V) 
Df dt 



(13) 



These are completed by the expressions for the radiation flux: 

H = -i^r 3 vr, (14) 



3i<p 



and for the specific entropy: 



(y-l)p \p?l 



(15) 



The physical quantities, which intervene in these equations, are 
the pressure P, the temperature T, the density p, the gravita- 
tional potential created by the star cD, the viscous force F v i SCl the 
nuclear energy production rate by unit mass e, the absorption 
coefficient k, the ratio y = Cp/Cy of specific heats at constant 
pressure (Cp) and constant density (C p ), the molecular mass 
p. cr is the Stefan-Boltzmann constant, and K g the perfect gas 
constant. 

Treating the tides as a small amplitude perturbation of the hy- 
drostatic structure of the object, we split all scalar quantities as 



X (r, 0, ip, t) = X (r) + X' (r, 0, if, t) ; 



(16) 



Xq designates the spherically symmetrical profile of X, and X' 
represents the perturbation due to the tidal potential. The refer- 
ence state is drawn from an up-to-date stellar or planetary struc- 
ture model. 

We now turn to the global velocity field W, from which we sub- 
tract the star's rotation (assumed here to be uniform) to isolate 
the time dependent tidal velocity field V : 

W(r,e,(P,/) = fiAr + V (r, 0, ip, t) . 



This tidal velocity field may be expanded in the same way as the 
tidal potential (6), and each component V , may be analyzed 
in the reference frame %c rotating with the corresponding tidal 
frequency, i.e. the angular velocity 



(17) 



In that reference frame, all quantities X' induced by the equilib- 
rium tide are stationary, in particular the density and the velocity 
field: 



Dp' 



Dl 



0, 



DV , 

m,l 



Re 



Dl 



0. 



(18) 
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From the continuity equation ( 10), we then deduce that 

V-V;„ >; = 0; (19) 

this property will hold in any reference frame, and thus for the 
total tidal velocity field V. This confirms the solenoidal char- 
acter of that flow, contrary to the claim made by Scharlemann 
(1981). 

3.2. Adiabatic system - Dissipative system 

In our treatment of the equilibrium tide, we shall separate all 
quantities and equations in two parts. 



First, an adiabatic system (I) in phase with the perturbing 
potential U (eq. 8) : it corresponds to the adiabatic tide, 
i.e. the star's response to the tidal excitation, ignoring all 
dissipative processes. 

Second, a dissipative system (II) in quadrature with the per- 
turbing potential : it corresponds to the star's response due to 
the dissipative processes (here, the turbulent friction due to 
the convective motions). 



stars and giant planets - The coplanar case 

We apply the same distinction to the adiabatic tidal flow: 

Y =Y l (r,e,<p,t) + V n (r,e,<p,t) , (23) 

where the tidal field Vi induces a radial displacement which is 
in phase with the potential U: 




Fig. 2. Typical system of tidal interaction. Body B exerts a tidal 
force on body A, which adjusts itself with a phase lag 5, because 
of internal friction. This adjustment may be split on an adiabatic 
component, which is in phase with the tide, and a dissipative 
one, which is in quadrature. 



The scalar perturbations (16) will thus be split as: 

X' (r, 9, ip, t) = Xj (r, 9, if, t) + X u (r, 9, if, t) , (20) 

where X\ designates the contribution of the adiabatic tide, in 
phase with the potential U : 

Xi = Y X u (r, 9) cos (2if - hit) 
i 

= P 2 (C ° S 0) COS (2 f ~ lni) ' ( 21) 

and X\\ that of the dissipative tide, in quadrature with U : 

X u = YjXuJ^^ sm(2if-lnt) 
i 

= ^ X nj (r) ?\ (cos 9) sin (2ip - Int) . (22) 



u\(r, 9) sin(2y; - Int) 
v\(r, 6) sin(2i,p - Int) 
w\(r, 9) cos{2if — Int) 



(24) 



and the velocity field Vn associated with the dissipative tide 
which is in quadrature: 



Vn = E 



un(r, 8) cos(2</; - Int) 
vu(r, 9) cos{2if — Int) 
- wu(r, 9) sin(2^ - Int) 



(25) 



We now linearize the original system (9 - 11), and split it like- 
wise in systems (I) and (II), both with their equations of motion, 
of continuity, Poisson, and of entropy, in this order. 



Adiabatic system (I) 
DV, 



Po 



Dl 



+ 2£2 A Vi 



-VP I -p V(<D I + C/)-p I V(Do, (26) 



Dp, 
Dt 



(27) 

(28) 
(29) 



Dissipative system (II) 

[DV n 



Po 



Dl 



+ 2il A V 



- VP„ - po VO„ - pa Vd> () + F u , (30) 



^+pV.Vn = 0, 
V 2 <£> u = 4^p u , 



Po- 



Dl 



= 0. 



(31) 

(32) 
(33) 



As was shown in Z66, the dissipative tide is much weaker than 
the adiabatic tide. That means that for a given scalar field X' = 
X\ + X\\, we have Xu <K X\. Moreover, this allows us to neglect 
the viscous force Fi(Vn) in system (I), but we keep Fn(Vi) m 
system (II), where it induces the dissipative tide. 

3.3. The adiabatic equilibrium tide 

3.3.1. Adiabatic response to the tidal potential - Love number 

We begin with the study of the adiabatic equilibrium tide, which 
is the hydrostatic response to the perturbing potential U exerted 
by the companion body, in the absence of dissipation. We expand 
the perturbations of pressure as in (21): 



Pi (r, 9, ip,i) = ) Pi,, (r) ¥ 2 2 (cos 9) cos {2xp - Int) , 



i 



(34) 
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and likewise for pi(r, 9, <p, t) and <E>i(r, 9, <p, t). These perturba- 
tions obey the equation of motion (26) of system (I), from which 
we remove the acceleration term on the l.h.s.: 



3.3.2. The poloidal component of the adiabatic velocity field 



VPt = -po V (fy + U)-pi VO 

From the vertical balance we draw 
dP d 

-j— = -Po rr ($1 + U)-pi go 
Or Or 

and likewise from the horizontal balance, 
Pi = -Po (<&i + tO , 



We have established earlier that each Fourier component of the 
tidal perturbation is stationary in the reference frame He rotating 
(35) w ^ tn ^ e an g u l ar velocity Q s = (1/2) n. This is true in particular 
for the potential (<t>i + U), and therefore 



(36) 



(37) 



D 
Dt 



(<Di + U) + Vi ■ VcD = 



(45) 



where go(r) = Qm(r)jr 2 designates the unperturbed gravity. 
Eliminating Pi we get the following expression for the density 
perturbation 



Pl = _l^ (a > I + [7); 

go dr 



(38) 



we implement it in the perturbed Poisson equation (28), taking 
into account that U is an harmonic potential (i.e. V 2 t/ = 0), to 
finally obtain (cf. Sweet 1950): 



From this we draw the vertical component of the velocity field 

Vi r = > o-, («, Q) U, P?(cos 9) sin(2co - Int) , (46) 

So V 

where its amplitude varies linearly with the tidal frequency 

o-, = In - 20 . (47) 

We have now to complete this vertical field with its horizontal 
counterpart. Since Vi is divergence free, with a radial compo- 
nent, it is natural to seek it in form of a poloidal field. Following 
Z66 and Rieutord (1987), we write 



v^(a, I + i/)-M^ ( o I + c/) = o 

go dr 



(39) 



Vf sVA VA[(r;r P )r| , 



(48) 



and expand the poloidal current function x P C*> <P> m lts 
The potential perturbation <t>i has the same dependence as U in Fourier components: 
(6, tp, t); therefore we may introduce a function h{r) which de- 
pends only on the radial coordinate, such that (<l>i + U) — hU. It 
obeys the following differential equation: 



X 



(x, 9, <p,t) = 'V f { (x) Pj (cos 9) sin (2tp - Int) 



(49) 



d 2 h 6 dh 

T1 + -- r +f(x)h = 0, 
ax L x ax 



(40) 



In the inertial frame the poloidal velocity field is thus given by 



where x = r/R, R being the total radius. The function fix) de- 
pends on the internal mass distribution at equilibrium: 



fix) 



4ngR dpo 
go dx 



(41) 



The solution of (40) is constrained by two boundaries condi- 
tions: regularity at center x — 0, and h — > 1 as x — > oo; it may be 
written as 

<E>i(jt, 9, <p,i) = V ®u(x) ?l (cos 9) cos (2ip - Int) 

i 

= [h(x)-l] U(x,6,<p,i) (42) 
= [h(x) - 1] V U, x 2 Pjicos 9) cos(2(p - Int) , 



6^fP 2 (cos0) sin (2tp -Int) 

1 d . 2 v dP 2 (cos0) 

xrt\ ) sin (2(p - Int) 

x dx v ' 

2 d 

x dx 



df) 



, „ „v P;(cos0) 

( x 3 r) i n c ° s (2^ - to) 



sinfl 



(50) 



Comparing the radial component with (46) we obtain the follow- 
ing expression for t$ (x): 



1 hx 2 

TO-, U,. 

6 go 



(51) 



where from (8) 



1 

Ui = ~- 



m B 



4 niA + nig 



n R G 2 ,o,i-2(e). 



(43) 



3.3.3. The toroidal component of the adiabatic velocity field 

The description of the adiabatic velocity field is not yet com- 
plete. Indeed, in the presence of rotation, the equation of motion 
(26) can not be satisfied by the sole poloidal velocity field, but 



The surface value of the ratio <£>\/U is the second-order Love one must add a toroidal field 
number in the adiabatic case: 



V[=V A ( x T r) , 



(52) 



(44) 



This constant is twice the classical apsidal motion constant, and 
it is approximatively the fifth moment of the mass distribution 
(see Eq. 2.15 of Zahn, 1966a). 



where x T is me toroidal stream function. Like x P before, this 
function may be expanded in Fourier series: 



X 



T (x, 9, if, t) = YjXl (x, 9) cos (2<p - Int) , (53) 
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with the following projections in spherical coordinates 



-^-.x] sin (2^ - hit) 
sin# ' 



" de 



cos (2<p - Int) 



(54) 



Here, our treatment departs from that followed in Z66. Instead of 
working in the inertial frame, we adopt the reference frame %c 
rotating with the Fourier component of the tidal potential poten- 
tial U, i.e. with the angular velocity Cl R — In. In this frame, all 
quantities are stationary; hence, the momentum equation reduces 
to 



po 2Q K A Vi = - VPi - po V (<Di + U)-pi V<D , (55) 
for each Fourier component. Their 6 and <p projections may be 



written 

2Q* cos 6 
and 



86 



P 2 (cos0) 2 d 



sin# 



x dx ^' 



2£2 



cos 
sin6> 



2x] + sin 6V\ (cos 0) 6ff 

dP 2 ( C OS<9) 1 d t , - 

cos 6^— -rtftf) 

d0 xdx 



dT(r,0) 
80 



2T(r,0) 
sin 6 



where the function <F includes all functions of the r.h.s. of that 
momentum equation. Its elimination, equivalent of taking the 
curl of (55), leads to 



2x <" 6 *'sr»d« 



— [sin 2 0P 2 (cosS)] 



, dP 2 
6 cos 9 P;(cos 0) + sin — - 
2 d<9 



(57) 



After some straightforward manipulations (see Z66), we get 
X T = X ?l W P 3 ( cos ) cos ( 2 V - • 



with 



d / hx 



-cr/x'-l- U,. 

15 dx\g 0l 



(58) 



(59) 



and reach finally the following expression for the toroidal field, 
in the inertial frame: 

T 2 i d / 



P 2 (cos 0) 

2—. 

sin0 

dP 2 (cos 6) 
d0 





sin (2<p - Ini) 
cos (2<p - Int) 



(60) 



Using the reference frame rotating with each component of the 
tidal potential, we are thus able to decouple rigourously the equi- 
librium tide, which is stationary in this frame, and the dynamical 
tide. Therefore, the 'pseudo-resonances' found in Z66 disappear. 



3.4. The dissipative equilibrium tide 

We proceed now with the treatment of the dissipation of the equi- 
librium tide. We assume that the scale separation (in time and 
size) between the large scale tidal flow and the turbulent con- 
vective motions is sufficient to allow their interaction to be rep- 
resented by an eddy viscosity j](r) acting on the tidal velocity 
field Vi. For simplicity, we assume this viscosity to be a scalar 
function of r, although we are aware that it is actually a diagonal 
tensor, as was shown by Penev at al. (2007). 1 This viscous force 
causes a lag in the tide and it projects it on what we have intro- 
duced above as system (II), giving rise to a flow Vn, which is 
in quadrature with the adiabatic tidal flow. For brevity, we shall 
call it dissipative velocity field. 

3.4.1 . The viscous force 

The treatment of dissipation of the equilibrium tide requires first 
to express the viscous force induced by the adiabatic velocity 
field Fn = F v i sc (Vi). Since Vi is solenoidal, that force takes the 
form: 



(56a) F n (V,) = V77 A (V A Vi) — V A V A (77VO + V (Vi ■ V77) - VjV 2 /? . 

(61) 

Once again we split the adiabatic velocity field in poloidal and 
toroidal components: 

F II (V I )=Fn(Vf)+F n (v I r ), (62) 

(56b) f° r which we find the following expressions in terms of the func- 
t; "™ fi f and t;J 



tions £f and we have introduced in (49) and (53): 



// (x) P 2 (cos 0) sin (2<p - Int) 



f?(x) 



dP 2 (cos 0) 



2f, H (x) 



d0 



P 2 (cos 6<) 
sin 9 



sin (2(p - Int) 



cos (2(p - Int) 



(63a) 



where 



f > R 2 \ x 2 



R 2 xf^n 



+2 



d7?d£f 
dx dx 



dj] 



(63b) 



dx 3 v ' ~ dx 



and 



dx 2 



2/, r (x) 



P 2 (cos 0) 
sin0 

dP 2 (cos 0) 
d0 



with 



f 



R 2 



sin (2ip - Int) 
cos (2(p - Int) 

( tT\ 



(64a) 



dr] d 
dx dx 



(64b) 



1 Other effects would have to be considered if the rotation were not 
uniform, as we assumed here; the case of differential rotation will be 
addressed in a forthcoming paper. 
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Fig. 3. Top-left: 3-D view of the total (poloidal and toroidal) adiabatic equilibrium tide velocity field (white arrows). The red and 
orange arrows indicate the direction of the primary's rotation axis and the line of centers respectively. Top-right: Representation of 
this velocity field at the surface of the primary (black arrows); the color-scaled background represents the normalized tidal potential 
intensity (blue and red for the minimum and maximum values respectively). Bottom-left: View of the velocity field (white arrows) 
in its equatorial plane of symmetry; the color-scaled background represents the velocity value (black and orange for the minimum 
and maximum values respectively). Bottom-right: View of the velocity field (white arrows) in its meridional plane of symmetry; 
the color-scaled background represents the velocity value (black and purple for the minimum and maximum values respectively). 



3.4.2. Toroidal velocities of the dissipative tide 

Let us rewrite here the momentum equation that governs the dis- 
sipative tide (cf. 30) in the reference frame rotating with the l- 
component of the tidal potential: 



Po 



[m R A V n ] = - VPn - po VO„ - pa V<D + F u . 



(65) 



The viscous force Fu has a toroidal component (64a) and it gen- 
erates therefore a toroidal velocity field Vj which obeys 



po[2n*AV£] = F n (vf) 



(66) 



To determine this velocity field, we introduce the stream func- 
tion D, (x, 0) and project on and <p: 







2D, (x, 0) cos (2<p - hit) 



-— [sin £>/(*, 0)] sin (2<z> - Z«f) 
o6 



(67) 



Next, we implement the Coriolis acceleration associated with 
this velocity field in (65), which we project on the horizontal 
coordinates: 



2p n R cos [sin GDi] = + 2 /, 



cT * 3 



(COS 0) 



do 



30 



sin 



4p cos 0D, = — + fj 

sin d0 



(68a) 



(68b) 



where T contains the scalar functions of the r.h.s. of the momen- 
tum equation other than the toroidal part of the viscous force. As 
before we proceed with the elimination of T and obtain this ex- 
pression for D,: 



D, = d,(x) 



(COS0) 



sin0 



, where d, = 



// r = 7— fl ■ (69) 
Inpo 



1 



F. Remus et al.: The equilibrium tide in stars and giant planets - The coplanar case 



Let us rewrite the toroidal velocity field induced by the dissipa- 
tion tide: 



Q*p 
3 



P 2 , (cos ( 



cos (2(p - Int) 



fl 



dP 2 (cos6») 
d9 



sin (2(p - Int) 



(70) 



The viscous force, as well as its associated toroidal velocity 
field, produces a mass redistribution inside the star or planet, 
and thus a perturbation of density pn and gravific potential <l>n- 
This causes a tidal torque to be applied to the object, which will 
be responsible for the dynamical evolution of the system. 

3.4.3. Perturbation of gravific potential 3> n 

We saw in the previous section that viscous dissipation gen- 
erates a toroidal velocity field Vj, which is in quadrature 
with the toroidal part of the adiabatic flow. It remains 
now to determine the gravific potential On, which is associ- 
ated with it, and which is in quadrature with the tidal potential U. 

Perturbed Poisson equation 

Let us go back to the Navier-Stokes equation (30) satisfied by 
the toroidal velocity field whose expression is given in (67), 
taking into acount relations (63a) and (64a) of the viscous forces 
due respectively to the poloidal and toroidal adiabatic velocity 
fields. Projection on r and 9 leads us to the following equations: 



2Q. R p {) sin 9— (sin 9 D,) sin (2ip - Int) = 
89 

1 dP u 1 3<D„ 1 cKDo fV u2 , ^ . „ , . 

-Po-^— -Pii-^— +/, P 2 (cos6>)sin(2^-/«f) , 
R ox R ox R ox 



(71a) 



and 



2O K p sin 9— (sin 9 D,) sin (2tp - Int ) = - - — — - — 

o9 R x 39 R x 39 



„dP 2 (cos60 T P 2 {cos9) 



sin# 



sin (2(p - Int) . (71b) 



The Coriolis acceleration on the l.h.s. of these equations projects 
on the associated Legendre polynomials of second and fourth 
order P 2 and P 2 A , since 



Q 

2Q*po sin6»— (sin 9 D,) 
o9 



d 

6 sine— P 2 (cos 6) ff 
^[-20P 2 + 6P 2 ]/, r . 



(72) 



Thus, the viscous force causes perturbations of density, pressure 
and internal potential that contain these spherical functions only, 
and we may expand them as 

X u = V [X 2 ,i (x) P\ (cos 9) + X 4J (x) P\ (cos 9)] sin (2tp - hit) . 

(73) 



However, we agreed earlier that the tidal potential could be rep- 
resented by a quadrupole only. In this approximation, the pertur- 
bations in P 2 do not contribute to the tidal torque, and we will no 
longer consider them. Keeping only the P^ harmonics, dropping 
the index II, and setting s — 2 , we rewrite (71) as 



~2Q. R p {) d, 

25 d 

— 2Q xp di 



1 dP, 
R~dx~ 



1 d(D, 
R dx 



P^poQt+Rxff 1 +5RxfJ 



(74a) 
(74b) 



where we have used the identity 2P 2 (cos#) = 5 sin#dP^/d#. 
By eliminating the pressure and gravity perturbations, we isolate 
the term, which will be implemented next in the Poisson equa- 
tion (32): 



gop, = -^o,+z,w , 

dx 



where Z, 



f 



d , „x 40 



dx 



which finally yields 
xd?^ 



-yZ 



AnQR 2 



(75) 



(76) 



where as above /(x) = -{AnQR/go) (dpo/dx). Itremains to solve 
that differential equation with the appropriate boundary condi- 
tions: 

d$, - - 

x + 3<3>/ = at x = 1, and <t>/ regular at x = 0. (77) 

dx 

Application to a thin convective envelope 

Equation (76) allows for an approximate solution in the case 
where the star, or planet, possesses a convective envelope. If that 
envelope is not too thick, we may assume that its mass is negli- 
gible, so that go(x) » g s /x 2 , and that the function h (cf. 40), may 
be assumed constant: h{x) as 1. Then the poloidal and toroidal 
stream functions take the simple monomial form 



1 



U, 



si = 7°"' — •* > 
6 g s 

>r 2 4 
Si =-z(Ti — x , 



(78) 
(79) 



where g s = g (l) designates the surface gravity. From (63b), we 
then deduce that: 



fV U ( rH\ U' 



dx 



and from (64b) 



R 2 gs 



: (24x 2 T] + 24x 3 rj' + 4x 4 i]") , (80a) 



dx 



R 2 g s 
2 



- [48x 2 r] + 29x 3 rj' + 3x 4 ?]"] . (80b) 



Then, introducing these expressions in (76), we get: 
Ui 



X - [936 x 2 r] + 632 x 3 r( + 76 x 4 77"] 



(81) 
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The last step consists in solving Poisson's equation (76), which 
we rewrite 

\jiQR 1 



d- 



dx 1 



3-7 

x Z/ 



-zi{x) 



(82) 



letting again go(x) ~ gs/x 2 and neglecting the function f(x) 
compared to 6/x 2 (for a justification of this approximation we 
refer to Z66). This can be done semi-analytically by using the 
method of Green's integral kernels (Morse & Feshbach 1953, 
see also Z66); the formal solution is 



<D,(x) 



(83) 



At the surface On,/, as we name it again to recall the dissipative 
origin of this tide, takes then the value 



<&„,,(!) 



AnQR 2 U, 

"357 
-l 



R 2 g, 



a- 1 



f [936 x 8 tj + 632 x 9 77' + 76 x 10 77"] dx , 



(84) 



where the summation is made over the convective envelope, 
whose lower boundary is at x = x^cz and the upper boundary 
at x — 1 . Thus the gravific potential due to viscous dissipation 
On/, which is in quadrature with the tidal potential £//, is pro- 
portional to the tidal frequency <x/, provided that the turbulent 
viscosity is not a function of cr/. 
Eq. (84) can be recast as 



%,/(!) , 2088 

= -47T ■ 



R- 



Ui 



35 m A g : 



f 



x S T](x) dx . (85) 



We see that this ratio, which quantifies the tidal dissipation, is 
a function of the tidal frequency, unless the turbulent viscosity 
is inversely proportional to that frequency. Finally, let us stress 
that the equilibrium tide dissipation depends strongly on the in- 
ternal structure of the considered body, and in particular on the 
thickness (or the mass) of the convective envelope. 

3.4.4. Turbulent viscosity 

As stated above, we assumed that the convective motions and 
the tidal flow are separated enough in temporal and spatial scales 
that their interaction can be described by an eddy viscosity 77 act- 
ing on the tidal velocity field. As long as this scale separation is 
fulfilled, we may take 77 = pov, (x) * po V t (x) I (x), where v, is 
the kinematic turbulent viscosity, V, the vertical rms velocity of 
the turbulent eddies, and I their mean-free path or mixing length. 
The radial profile of 77 is readily drawn from a classical mixing- 
length model, since jf = (3/20) a A p\ Hp (x) 3 F com (x), where 
Hp = |dlnPo/dlnr| _1 is the pressure scale-height, a = £/Hp 
the mixing-length parameter and F conv the convective enthalpy 
flux (cf. Z66, Zahn 1989). 

However, this simple prescription no longer holds when the 
tidal period P^ = 2n/cri becomes shorter than the life span 
fconv = (/Vt of the convective elements. One expects then the 
tidal dissipation to be reduced: just how much is still a matter of 
debate. On phenomenological reasons, it seems then natural to 
replace the mean-free path by the distance traveled by the con- 
vective eddies during, say, half a tidal period; then the eddy- 
viscosity varies proportionally to the tidal period (Z66): 



v t = \v}Pmz =nVf(Tj l . 



(86) 



Recently, this scaling has been confirmed through numerical 
simulations performed by Penev et al. (2007) with a code de- 
signed to model stellar convection. They found also that v, is in 
fact a diagonal tensor, with stronger transport in the vertical than 
in the horizontal direction, but we shall not take into acount these 
refinements and simply treat 77 as a scalar. 
Another prescription has been proposed by Goldreich & Keeley 
(1977), where the turbulent viscosity varies as the square of the 
tidal period. It implies a much stronger reduction of the viscos- 
ity that is incompatible with most observations, for example the 
circularization of binaries with a red giant component (Verbunt 
& Phinney 1995); but it could be adequate for the dynamical 
tide involving inertial modes, which develop boundary layers 
that are of much shorter scale than that of the equilibrium tide 
(cf. Ogilvie & Lin 2004). 

It thus appears that turbulent dissipation operates in two regimes, 
depending on how the tidal period compares with the local con- 
vective turn-over time, which in a convection zone can vary with 
depth by several orders of magnitude. To ensure a smooth tran- 
sition between these two regimes, one may take 



v t = v t e 



1 + 



2/c 



-1/2 



1 + 



(^conv V 
n I 



-1/2 



(87) 



as illustrated in Fig. 4. 

In the upper part of a convective envelope, where the convective 
turnover time is shorter than the tidal period, v, does not depend 
on the tidal period; the tidal dissipation varies proportionally to 
the tidal frequency (cf. Eq. 85). But in the opposite case, when 
the life span of the convective eddies exceeds the tidal period, 
which is likely to occur at the base of convection zones, tidal 
dissipation is independent of the tidal frequency. Note that these 
two regimes still persist once the summation over 77 has been 
performed in Eq. (85). 



4. The dynamical evolution of binary systems 

Now that the tidal velocity field and the associated mass redis- 
tribution have been obtained, we are ready to determine the tidal 
torque, which causes the exchange of angular momentum be- 
tween each component and the orbital motion. Let us recall the 
role played respectively by the adiabatic and the dissipative tide. 
First, the companion exerts a tidal force on the considered com- 
ponent, which is elongated in the direction of the line of centers: 
this mass redistribution described by the perturbed density p\ 
corresponds to the adiabatic tide. In the absence of dissipative 
processes, the induced torque Fi = dU/dippi dfV presents pe- 
riodic variations of zero average and no secular exchanges of an- 
gular momentum are allowed. However, turbulent friction gener- 
ates a mass redistribution described by pn. This dissipative tide 
produces the torque Fn = dU /dippn d'V of non-zero average, 
which causes the net exchange of angular momentum between 
the components of the system as illustrated in Fig. 2. 

4.1. The dissipation quality factor 

The expression of the torque Fn involves the ratio On,//£// which 
we have established above in Eq. (85), or equivalently the quality 
factor Q which we shall now introduce. For this it is convenient 
to work in complex notations, with the real part representing the 
adiabatic tide and the imaginary part characterizing the dissipa- 
tive contribution. For the perturbed gravific potential we have 
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then 



<£' ( r , t ) = %: {o ? (x)P|(cos 6) exp [i (2<p - bit)] } , (88) 

where we define, following the notations of Eqs. (21-22), 

O, ^Oy-iOn,,. (89) 

For each Fourier component I, the complex second-order Love 
number k 2 , which characterizes the deformation response of A 
to the stress exerted by B, is defined by 



k 2 (cr,) = 



OKI) 
Ui 



(90) 



Since the turbulent friction introduces a phase delay ee [ _ f > f ] 
of the reaction relative to the load, k\ is also of the form 



h (cr/) = k 2 (cr/) e 



-lV(cr ; ) 



(91) 



where k 2 is the norm of k 2 . From Eq. (90), we identify the real 
part of k 2 as the adiabatic second-order Love number given by 
Eq. (44), thus 



kf 



<&(k 2 ) 



1>y(l) 
Ui 



k 2 cose. 



(92) 



The ratio <t>n,;(l)/t// = k 2 sin e is similar to k* d = k 2 cos e, but ap- 
plied here to the dissipative tide. By analogy with an electric cir- 
cuit, we then define the quality factor Q > (see Ferraz-Mello 
et al. 2008; Efroimsky 2012) by 



1 



\lm k 2 (cr/)| 



Qicrd k 2 (cr,) 
such that the tidal dissipation is given by 



= sin \e(cri) \, 



(93) 



k 2 (<ri) 
Q(cr,) 



= \lm k 2 (cr/) = 



®ny(l) 



An- 



2088 R 2 
35 m A g s 



Ui 



f' 



x i](x) dx 



(94) 



where we have recall the expression of Eq. (85). Note that the ra- 
tio k 2 (<Ti)/Q (cr/) can be frequency dependent, whereas k^ d does 
not. This tidal quality factor is also linked to the geometrical lag 
angle 5, which is half the phase lag e (see Efroimsky & Lainey, 
2007): 



k 2 ((Ti) 



= k 2 (cri) sin |e (cr,) | = k 2 (cr ; ) sin \26 (cr,)\ 



(95) 



Recalling Eq. (85) of <t>n,/(l)/t/,, we notice that when the turbu- 
lent viscosity does not depend on tidal frequency, the tidal lag 
angle 6 (cr/) = cr/ Af (cr/) is proportional to the tidal frequency, 
which means that the time lag Af(cr/) of the tide is then con- 
stant, and takes the same value for each Fourier component of 
the tide: this is what has been called the weak friction approxi- 
mation (Hut 1981). It corresponds to the slow tide (? conv < Aide), 
as illustrated in Fig. 4. In the opposite case where f conv > Aide, 
the tidal angle 8 is independent of the tidal frequency, and so are 
also the phase lag e and the quality factor Q, such that the time 
lag Af = (5/cr, is inversely proportional to the tidal frequency. 




Fig. 4. The two regimes of turbulent dissipation (Eq. 87). As 
long as the local convective turn-over time remains shorter than 
the tidal period (t com < Aide), the turbulent viscosity v, (in black 
dashed line) is independent of the tidal frequency, and the inverse 
quality factor k 2 /Q = \<S>n/U\ (in red continuous line), given by 
Eq. (94), varies proportionally to the tidal frequency cr/ (so does 
also the tidal lag angle). When f conv > Aide, v t varies proportion- 
ally to the tidal period, whereas k 2 /Q — \<f>u/U\ does no longer 
depend on the tidal frequency. Note that v, and k 2 /Q have been 
scaled by the value they take respectively for f CO nv/Aide — * and 

+oo. 



4.2. The disturbing function 

Here, we shall use the disturbing function method as described 
in Brouwer & Clemence (1961), Kaula (1962), Yoder (1995-97), 
Correia & Laskar (2003a,b) and Mathis & Le Poncin-Lafitte 
(2009) (hereafter MLP09) to determine Fn and the associated 
variation of the bodies' angular momentum and of the orbital 
keplerian elements of B (namely, a and e). 

The first step to achieve this goal is to introduce a supplemen- 
tary body C with a mass mc, which we call the orbiter, to the 
primary A and the perturber B (see Fig. 5 and Kaula 1962). The 
keplerian elements of its orbit are a*, e*, and M* = n*t with «* 
the associated mean motion. 

Next, we identify the disturbing function on C to the perturbed 
external gravific potential of A: 

% = <D' (r c , i) = Yj* i® 1 ' (xc) P 2 (0) exp P (2 ^ c ~ fw * f)] } ' 

(96) 

where rc = (rc =xcR,6q =7t/2,i/>c), with xq > 1, is the 
position of C in the equatorial plane of A, which is also the 
plane of the orbit, and where the index I* is associated to the 
Fourier expansion of the orbital motion of C. 

From now on, we choose to follow the systematic method de- 
veloped in MLP09 to derive the expression of % as a function of 
the keplerian elements of B and C. In this work devoted to the 
study of the tidal dynamics of extended bodies, the perturbing 
function is derived for components X having an external gravific 
potential written as 



M, 



J X =0 Hl X =-iX 



'x 



(97) 
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expression, once applied to our case: 



Q 



iCF,.->P,)| 



Fig. 5. Here, we illustrate the used method of the perturbing 
function. The primary A is deformed because of the tides exerted 
by the perturberB. Because of A's internal friction processes, the 
tidal bulge is shifted from the line of centers with the tidal angle 
Si given in Eq. (94) for a given I Fourier component of the tidal 
potential. The keplerian elements of B's orbit are a, e, and M. 
Next, a third body C, which orbits around A with a mean motion 
n* , is introduced, and the variations of the keplerian elements of 
its orbit (a*,e*,M*) are derived. After using Lagrange's equa- 
tion (Eqs. 108-109), we finally consider that the orbiter and the 
perturber are the same (i.e. C=B). 



where (rx, Ox, <Px) are the spherical coordinates attached to their 
equatorial planes. We introduce the gravitational moments of X 



M.„ 



An C . 
2s X + 1 Jy x ' 



~sx ym x 
'x 1 s x 



(0x,7px)pxdV x , (98) 



where dVx = ^d?x sin 6x dOxdipx', the non-spherical compo- 
nents of O x are due to tides and to the own structural shape of 
the body. We use spherical harmonics 



A / g£±l^ ■ gV (cose)exp[imy] 



\ An (s + \m\)\' 



(99) 

Here, we apply this general formalism to the interaction between 
A, which has a tidal gravitational moment M^ 2 (the quadmpolar 
approximation is assumed for the perturbed external gravific po- 
tential of A as for the tidal potential U), and C, which is assumed 
to be ponctual as well as B. Then, the Kaula's expansion (Kaula 
1962) allows to convert the expansion on spherical harmonics on 
to a Fourier one on keplerian elements (see Eq. 79 in MLP09). 
Applied to our coplanar case, this gives for X = {B, C) 



Y^n/2, Vx ) 



where 




2 G 2,o,; x -2(ex)exp[i>l', x ] , (100) 



% x =IxMx + 2ojx-2®a 



(101) 



The spherical coordinates of the X position in the A refer- 
ence frame are (rx,7i/2,<px) while lx = {1,1*}, ax = {a, a*}, 
e x = \e, e*), and M x = n x t = \M - nt, M* = n*t). The argument 
of the periapsis lox = {0, of) is also introduced since it could 
be different for B and C. Finally, ©a is the mean sideral angle 
between the minimal axis of inertia and the (AX) axis, and thus 
Q = d® A /df. 

This transformation leads in the general case to the expansion of 
% given in Eq. (144) in MLP09, which reduces to the following 



(102) 

where X,i (o~i) is the so-called dissipative impedance, defined as 

~ ( , him) , $n./(l) , .. , . 
Xi(cri) = —. =1-1— = 1—1 tan e (<x ;) , (103) 



while 



15 



1/2 



H x (ex) = | — G 2A , x -2 (ex) ■ (104) 



i(T,.-f,) I 



Thus, we can recast Eq. (102) as 

*= - j> ih (m) 4 i^S-n WW e^ 

I p \ J a a 

(105) 

The adiabatic component Hi (which only involves "Re (£2)) of 
the perturbing function % contributes only as zero-average per- 
turbations to the angular momentum and the orbital elements, 
whereas its dissipative component 'Rw (i.e. the one proportional 
to Im (k-2p given by 



ki(o-i) An Qm B R 5 

ft\Q((Tl) 5 fl 3 fl *3 



n(e)n-(e*)sm(y,, - V,) 



v [ Oiu(l) An Qm B R 5 fVI w ,\ 

= -h\—Tl^ n{e)n ~ (e )sin( ^ - %) ] 

(106) 

acts on secular time-scales. 

We thus obtain for the mean disturbing function the same form 
as Correia & Laskar (2003a). The main difference is that in our 
case we link the viscous dissipation with the physical processes 
operating in the considered bodies,i.e. here turbulent convection. 
Since the turbulent viscosity depends on the tidal frequency, we 
have to treat the dissipation separately for each mode / of the 
tidal potential, as it has been already shown in Zahn (1977) and 
Ogilvie & Lin (2007). The general method presented here, based 
on Zahn (1977) and MLP09, thus allows to go beyond the usual 
constant Q or tidal lag approximations that can only constitute a 
first approach as shown in §3.4.4 and in Ogilvie & Lin (2004). 
In a near future, we will be able to apply this method to the fluid 
dynamical tide and to the anelastic tide for planets that have a 
rocky part (Efroimsky & Lainey 2007; Henning et al. 2009). 

4.3. Equations of dynamic evolution of the system 

We are now ready to derive the evolution equations for the an- 
gular momentum of A and for the orbital elements of C. First, 
following Yoder (1995-1997) and Correia & Laskar (2003a,b), 
we obtain the equation for the angular momentum of the primary 



df 



dI A Q. 
df 



m c d-Q %i, 



(107) 



where I A designates its inertial momentum. Next, applying 
Lagrange's equations (see for example Brouwer & Clemence 
1961), we get those for the evolution of the semi-major axis of 
the C's orbit and its eccentricity: 



da* 2 „ „ 

-77 = —dwRn , 
df n a 



(108) 
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de* Vl -e* 2 1 - e* 2 

at n a re n a re 

Finally, variations of its mean anomaly M* given by 

dM* 2 1 - e* 2 

n* - —d a -<Rn - - . de.Hu (1 10) 



dt 



n* a* 



n* a* 2 e* 



are filtered. 



Applying previous Eqs. (108-109) to Ku (Eq. 106) using results 
derived in MLP09 (Eqs. 134-136-137), we get assuming that the 
orbiter is the perturber (i.e. B=C) 



dL A %n&m 2 B R 5 ^ | <D U ,/(1) 



dt 



5 a 6 



z 



Ui 



[H(e)f 



(111) 



where we identify Tn the secular torque applied to A, and 
1 da 2 4n@m B R 5 ^ L*n,/(1) 



a dt n 5 a 8 



S' 



Ui 



[Hie)]' 



(112) 



1 de _ 1 1 - e 2 4n Qm B R 5 
e dt n 



2 1 



5 a 
1 



+ (/-2) 



where we recall that (Eqs. 85-94): 



(113) 



%,/(!) 



= 4;r- 



2088 R 2 

35 OTAg, 



*J Art 



*BCZ 

(114) 

From these equations one readily draws the characteristic syn- 
chronisation and circularisation times: 



1 



1 



1 dL A 



(Q.-n)I A dt 



and 



1 



1 de 
e dt ' 



(115) 



Let us now comment this system of equations. First, as it has 
been already discussed in the previous section, the equilibrium 
tide dissipation is consistently described with its dependence 
on the tidal frequency. Moreover, we never made any assump- 
tion concerning the value of the eccentricity. More precisely, we 
chose to keep the most general expansion in e of the tidal po- 
tential U and the resulting disturbing function Ku. This means 
that, provided we take a large enough number of modes in the 
tidal potential expansion, as illustrated in Savonije (2008), we 
are able to treat the general case of eccentric orbits. Therefore, 
the present formalism allows to take into account the complex- 
ity coming both from a general eccentric orbital motion and from 
the physical processes that act on the tidal velocity field it gen- 
erates. 

Finally, let us remark that the equilibrium states resulting from 
Eqs. (111-112-113) will differ from those predicted assuming a 
constant Q factor or tidal lag because of the different dissipation 
of each Fourier component of the velocity field; this will be ex- 
plored in paper II where inclination and obliquity will be taken 
into account. 



5. Conclusion 

The originality of this work is the formulation of the equilibrium 
tide problem in the rotating frame associated with each Fourier 
component of the tidal potential. This allows to filter correctly 
the dynamical tide (due to the free oscillations of the considered 
body), and to avoid the so-called 'pseudo-resonances' which ap- 
peared in Z66 when the equations were written in an inertial 
frame. Morever, the solenoidal behavior of the velocity field 
of the equilibrium tide, which was questioned by Scharlemann 
(1981), is confirmed. We revisit the derivation of the adiabatic 
and the dissipative components of the equilibrium tide (the ve- 
locity fields and the associated perturbations of the gravific po- 
tential) with taking into account the most general expansion of 
the tidal potential in the case where the orbital spin and those of 
each body are aligned. This allows to treat the general eccentric 
case by taking a sufficient number of Fourier components into 
account. We obtain the associated rate of dissipation of the tidal 
kinetic energy into heat, given by the k^lQ ratio, which depends 
on the thickness of the convective envelope and on the tidal fre- 
quency (<t/). Indeed, because of the dependance of the turbulent 
friction on <x/, two dissipation regimes are identified. First, as 
long as the local convective turn-over time remains shorter than 
the tidal period (f conv < Aide), the turbulent viscosity v, is inde- 
pendent of 07, and the inverse quality factor k^/Q varies propor- 
tionally to the tidal frequency (so does also the tidal lag angle). 
When fern, > Aide, v, varies proportionally to the tidal period 
(cc <x~l), whereas kilQ does no longer depend on the tidal fre- 
quency. This emphasizes the need of going beyond the classical 
constant tidal lag approximation, not only for the dynamical tide, 
as it was shown Ogilvie & Lin (2004), but also for the equilib- 
rium tide. 

In the forthcoming paper, we shall take into account the orbital 
inclination and the obliquities of the bodies. Moreover, we will 
examine the impact of differential rotation on the equilibrium 
tide, and that of an anelastic core, which may well be present 
in giant planets. We plan also to check theoretical predictions 
againt observations in our Solar system as well as in exoplane- 
tary systems and binary stars. 
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